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We report on our recent results for deep-inelastic neutrino(V)-proton(P) scattering. 
We have computed the perturbative QCD corrections to three loops for the charged 
current structure functions F2, Fl and F3 for the combination vP — vP . In leading 
twist approximation we have calculated the first six odd-integer Mellin moments in 
the case of F2 and Fl and the first six even-integer moments in the case of F3. As a 
new result we have obtained the coefficient functions to 0(o? s ) and we have found the 
corresponding anomalous dimensions to agree with known results in the literature. 

1 Introduction 

In our recent research [1], we extended the program of calculating higher order perturba- 
tive QCD corrections to the structure functions of charged current deep-inelastic scatter- 
ing (DIS). Our studies are motivated by the increasingly accurate measurements of neu- 
tral and charged current cross sections at HERA with a polarized beam of electrons and 
positrons [2-4]. At the same time we are also able to quantitatively improve predictions 
for physics at the front-end of a neutrino- factory, see e.g. Rcf. [5]. To be specific, we con- 
sider neutrino(z^)-proton(P) scattering in the combination vP — vP 1 which corresponds to 
charged lepton-proton DIS as far as QCD corrections are concerned. Following Refs. [6-10] 
we compute the perturbative QCD (pQCD) predictions to three-loop accuracy for a number 
of fixed Mellin moments of the structure functions F 2 , F L and F 3 . 

Within the framework of the operator product expansion (OPE), and working in Mellin 
space, F2 P ~ vP and p^ p ^ vP are functions of odd Mellin moments only, while only even 
moments contribute to F^ p ~ vP . This is opposite to the case of the neutral current structure 
functions where only even Mellin moments contribute, and to the charged current case for 
vP + vP scattering [11], which is defined through the OPE for odd Mellin moments only. 
In the latter results for F^ p+vP and F^ p+vP to three-loops can also be directly checked in 
electromagnetic DIS [12,13] while parameterizations for F^ p+uP to three-loop accuracy are 
given in Ref. [14]. 

2 General formalism 

We consider unpolarized inclusive deep-inelastic lcpton-nucleon scattering, 

l{k) + nucl(p) -> i'(fc') + X, (1) 

where l{k), l'(k') are leptons of momenta k and k', nucl(p) denotes a nucleon of momen- 
tum p and X stands for all hadronic states allowed by quantum number conservation. In 
our research we are concentrating on charged current neutrino(^)-proton(P) scattering, i.e. 
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vP, vP via boson exchange. As is well known, the differential cross section for the 
reaction (TTJ) can be written as a product of leptonic L^„ and hadronic tensors 

da oc L^Wpv , (2) 

with L^ v for electroweak or pure electromagnetic gauge boson exchange given in the litera- 
ture, see e.g. Ref. [15]. The hadronic tensor W^ v in Eq. @ can be written in terms of so 
called structure functions Fi, i = 2, 3, L. 

We are interested in the Mellin moments of structure functions, defined as 



F{n,Q 2 ) 



dxx n ~ 2 F,(x,Q 2 ) , i = 2,L 



(3) 



and for F^{n, Q 2 ) one has similar relation with n replaced by n + 1 on the r.h.s. of Eq. ([3]). 
Here Q 2 = —q 2 > 0, q = k — k' and x is the Bjorken scaling variable defined as x = Q 2 / (2p-q) 
with < x < 1. 

With the help of the optical theorem and Cauchy's theorem from complex analysis one 
can relate the Mellin moments of structure functions to the parameters of the OPE for the 
nuclcon forward Compton amplitude T^ v : 



and the OPE for T„ v reads as 



i = 2,3, L 



(4) 
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^nuci ( n > A* 2 ) + higher twists , 



(5) 



where higher twist contributions are omitted. Ci jns denote the Wilson coefficients which are 
calculable in pQCD and A^ nci are matrix elements of quark non-singlet operators. The latter 
are not calculable in pQCD, rather they have to be extracted from experimental data. We 
restrict ourselves to quark non-singlet (ns) operators only since only these give nonvanishing 
contributions in the combination vP — vN (see Ref. [11] for details). 

Eq. ((4]) provides the basis to obtain Mellin moments of DIS structure functions in our 
approach by means of the OPE and the optical theorem. Furthermore, from the careful 
examination of the symmetry properties of the forward Compton amplitude and, related, 
the underlying Feynman diagrams, one can convince oneself that for the charged current 
vP — vN DIS, one encounters functions of only odd n for F2 and Fl and, functions of only 
even n for F3, respectively [11]. 

The pQCD calculation of Wilson coefficients Ci jns proceeds through the following steps. 
From the first principles we calculate the partonic forward Compton amplitude t„ v . The 
partonic equivalent of the OPE Eq. ([5]) for contains the same coefficients Ci !ns as in 
Eq. ([5]) and quark matrix elements ^4q S . Projection on the n'th Mellin moment of OPE and 



on the i'th parton invariant (i = 2, 3, L) with the help of the operator V^' i we get 
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Both sides of Eq. ^ are renormalized. In particular the renormalization of the local quark 
operator matrix element A™ gives rise to the factor Z ns on the r.h.s. of Eq. ([B]). This 
equation is our starting point for an iterative determination of the coefficient functions 
Ci.n S and the anomalous dimension 7 ns . The latter appears in a series expansion of Z ns in 
powers of the strong coupling a s and negative powers of the parameter e of dimensional 
regularization, D = 4 — 2e. The Ci !ns on the other hand are expanded in a s and in positive 
powers of e. Thus the l.h.s. of Eq. ^ leads to a well defined determination of Ci jns and Z ns 
in pQCD. 

3 Calculation and checks 

In the previous section, we have briefly explained the method to obtain Mcllin moments 
of the DIS charged current structure functions F^ p ~ vp , F^ p ~ pp and F^ p ~" p together 
with their respective coefficient functions and anomalous dimensions. To that end we have 
calculated the Lorentz invariants of the parton Compton amplitude t j ns , i = 2, 3, L, as given 
in the l.h.s. of Eq. ©. Due to the large number of diagrams involved in the calculations 
up to order a^. sufficient automatization is necessary. First of all, we have generated 3633 
diagrams up to three loops with the program Qgraf [16]. For all further calculations we 
have relied on the latest version of the symbolic manipulation program Form [17, 18]. 

For the treatment of Qgraf output, such as analysis of the topologies, the explicit im- 
plementation of Feynman rules etc. we have adapted a dedicated Form procedure conv.prc 
from previous work, e.g. Ref. [13]. Most importantly, this procedure tries to exploit as many 
symmetry properties of the original Feynman diagrams as possible in order to reduce their 
total number. 

For the calculation of the color factors for each Feynman diagram we have used the 
Form package color. h [19]. The actual calculation of the Mellin moments of the Feynman 
integrals has made use of the FORM version of Mincer [20]. Finally, on top of Mincer 
and Minos [7] some shell scripts managed the automatic runs of both programs for different 
parts of the calculation. 

We have performed various checks on our computation. Most prominently, we have kept 
all powers of the gauge parameter £ throughout the entire calculation for Mellin moments 
n < 10 to check that any ^-dependence vanishes in our final results. The Mellin moments 
with n > 10 were calculated without gauge parameter to facilitate the computations which 
become increasingly more complicated for higher Mellin n values. For these moments we 
also used TFORM [21], the multi-threaded version of Form. On machines with multi-core 
processors this leads to a significant speed up of our calculations, e.g. a speed-up of ~ 5 on 
a two-core four processor machine. 

We agree with the literature as far as the two-loop coefficient functions [22-26] and the 
three-loop anomalous dimensions [27] are concerned. In addition, for the first Mcllin moment 
of the coefficient function C2.11S wc have obtained exactly C2,ns — 1 to all orders in a s which 
is in agreement with the Adler sum rule for DIS structure functions, 

1 

J ^(f^{x,Q^)-F^{x,Q^=2. (7) 


The Adler sum rule measures the isospin of the nucleon in the quark-parton model and 
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does not receive any perturbative or non-perturbative corrections in QCD, see e.g. Ref. [28]. 
Therefore, this result is another important check of the correctness of our results. 

4 Conclusions 

We have reported on new results for Mellin moments of the charged current DIS structure 
functions F^ p ~ vp ', p^ p ~ vP an( j p" p ~ vP including the perturbative QCD corrections to 
three loops. In the former case (F2, Fl) we have computed the first six odd-integer Mellin 
moments while in the latter case (F3), the first six even- integer moments have been obtained. 
The results for F^ p ~ uP n = 1,3,5,7,9 and for Fg p ~ pp n = 2,4,6,8,10 are available in 
Ref. [11]. Results for n = 11 in the former case and for n = 12 in the latter will be published 
elsewhere. Finally, the discussion of phcnomenological consequences of our Mellin space 
results along with approximate parameterizations the coefficient functions in x are deferred 
to Ref. [29]. 
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